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Abstract 

In three-dimensional systems of the Ising universahty class the ratio of corre- 
lation length amplitudes for the high- and low-temperature phases is a universal 
quantity. Its field theoretic determination apart from the e-expansion represents a 
gap in the existing literature. In this article we present a method, which allows 
to calculate this ratio by renormalized perturbation theory in the phases with un- 
broken and broken symmetry of a one-component (/>^-theory in fixed dimensions 
D = 3. The results can be expressed as power series in the renormalized coupling 
constant of either of the two phases, and with the knowledge of their fixed point 
values numerical estimates are obtainable. These are given for the case of a two-loop 
calculation. 



1 Introduction 

The application of field theoretic methods to critical phenomena and second order phase 
transitions is one of the main tools for their quantitative theoretical investigation. Based 
on general theoretical arguments it is believed that physical systems, which undergo a 
second order phase transition, fall into distinct universality classes characterized by the 
dimensionality D of space, the number n of components of the order parameter and the 
underlying symmetry group. This manifests itself in the fact that within these universality 
classes many interesting quantities such as critical exponents and amplitude ratios have 
the same values |[|, |^ . 
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As a consequence of the long range nature of the effective interaction (for example 
between the spin variables of a ferromagnetic system showing a collective behaviour in the 
vicinity of the critical temperature T^) the microscopic structure of the given system can 
be neglected by introducing a real valued order parameter field = (0i(x), . . . , 0„(x)), 
X G R^, which is assumed to describe the phase transition. In this way one is led to 
a field theory with 0^-self- interaction with spontaneous breaking of its 0(n)-symmetry, 
which is one of the standard topics in quantum field theory. 

In the past many perturbative calculations have been made to estimate the values 
for the above mentioned critical quantities. To this end the e-expansion, initiated by 
M.E. Fisher and K.G. Wilson and subsequently elaborated by E. Brezin, J.C Le 
Guillou, J. Zinn- Justin and others for the massless 0^-theory in D = 4 — e dimensions, 
has been applied and yielded very good results. 

A quantity whose numerical value is important in comparing various experimentally 
determined universal quantities with theoretical predictions (see e.g. 0) is the univer- 
sal ratio of correlation length amplitudes /+//_. It is defined by the behaviour of the 
correlation length ,^ as a function of the temperature T near T = through 

. f f+t"" , t>0 T-Tc 

^ ~ < ; t:= (1) 

where the symbol ~ denotes the critical behaviour as defined in the usual way [|T], 0. 
The e-expansion of /+ / /_ is known up to second order for the second-moment correlation 
length and up to first order for the 'true' correlation length ^, which limits the accuracy 
of numerical estimates. 

In this article we follow an idea owing to G. Parisi 0, who suggested the use of 
renormalized massive perturbation theory in fixed dimensions D = 3. For the symmetric 
phase there exist extensive calculations by G.A. Baker, B.G. Nickel, D.I. Meiron et al. [|, |^, 
p!0| , who give the renormalization group functions and the critical exponents respectively 
up to six- and seven-loop order. Later C. Bagnuls, C. Bervillier et al. JIT], have 
developed a renormalization scheme, which among others makes it possible to determine 
the amplitude ratios of the susceptibility and the specific heat from these series. In general 
the quality of results from this method appears to be better than from the e-expansion. 

Up to now an extension of this method to the universal amplitude ratio /+//- is 
not found in the literature since this requires explicit calculations in the phase of broken 
symmetry, too. 

In order to fill this gap we consider Euclidean 0'^-theory with D = 3 and n = 1, which 
is a super-renormalizable field theory, that is believed to lie in the same universality class 
as the three-dimensional Ising model. Its Lagrangian density in the symmetric phase 
(t > 0) is given by 

CiM = ^(a</.o+(x)f + V(</.o+) 
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V(0o+) = lrnl+<Pl^{x) + ^go4A^) (2) 

with go > and mg+ > 0. The quantities in different phases are distinguished by indices 
+ and — , if necessary. In the broken phase {t < 0) the Z2-symmetry (p —0 is lost for 
the Lagrangian 

C{<p) = + V(0) 

V(0) = -^m^-0^(x) + ^^o0^(x) + ^^ = i, 0(0^(0:) -.of (3) 



with rriQ- > and the classical potential minima at = ±fo := ±y3mQ_/(7o- After an 
expansion of the potential around the positive minimum, which is equivalent to a shift in 
the field variable via 

(f)o^{x) := (j){x) -vo, (4) 

we arrive at 

1 1 / — 1 

V(0o-) = ^rnl-(j)l-{x) + —y3gomo^^l-{x) + —go^o-{x) . (5) 

These Lagrangian densities are the starting points of our perturbative calculations in 
both phases. It should be noted that the additional 0^-self-interaction of the field 0o- in 
the broken phase gives rise to tadpole-diagrams as well as to a non-vanishing one-point 
function, which is the vacuum expectation value of the field 00- . 

The Feynman rules in momentum space follow from the Lagrangians (§) and (|^) as 
usual. In particular each four- vertex is associated with a factor —go, each three- vertex 
gets a factor — i/S^mo- and the Feynman propagator is 

The divergences of the theory, appearing in the form of the only two primitively divergent 
graphs of the two-point function, are isolated by dimensional regularization. 



2 Renormalized perturbation theory 



Before introducing the renormalization schemes for the phases with unbroken and broken 
symmetry let us briefly outline our further strategy. 

To begin with we define the correlation length for general D as the second moment 



jd^xx'Gf^'\x) a^G(2.o)(p) 



(2,0), 



X) 



(7) 



of the connected two-point correlation function 

:= {Mx)My)) - (0o(x))(0o(z/)) . (8) 

Because of its relation 

-rrw = (G^°Hi>))" (9) 

to the two-point vertex function Fg^''''' in momentum space one verifies the following 
identity between the correlation length and the renormalized mass m^j: 



(10) 



As t'q''^\p) is given by 



-T^^''\p) = A-\p)-^{p), (11) 

where S(j9) is the sum of all one-particle irreducible two-point graphs with amputated 
external legs, equation (ITO) shows how C, can be determined perturbatively. Alternatively 



one defines the 'true' correlation length by the exponential decay of the two-point function. 
Then ^ equals the inverse physical mass (see e.g. ||I3|). 



From renormalization group theory the critical behaviour of the model is known to 
be controlled by the dimensionless renormalized coupling constant ur and its non-trivial, 
infrared-stable fixed point u*r. The calculation of renormalized quantities sketched in this 
section enables us to find analytic functions in upt, in terms of which the desired ratio 
/+ / /_ can be expressed, as will be shown in the next section. 

The renormalization scheme is established in terms of the bare (n, /)-vertex functions 
rS)"''H{P; "^o, go), with {p; q} = {pi, . . . gi, . . . , gj. They emerge from the expec- 
tation values (^(j)o{xi) ■ ■ -(poixn) |0o(?/i) ' ' ' |0o(l/O) after Legendre and Fourier transfor- 
mation. The renormalized mass (|10D is written as 



2_ (2,0). . 1 ._ dT^o'''\p-mo,go) 
T^R — -^3i-o lU,mo,5'oj, ^ 



(12) 



and in addition we define 



1. := -T^^''\{0-oy,mo,go) = -^T^^''\0;mo, go) ■ (13) 
Z2 orriQ 

The renormalization constants and Z2 play a prominent role in the theory because 
they stand in close connection to renormalization group functions and critical exponents. 
For the renormalized coupling we use different settings in each phase. 



Symmetric phase 

As usual |1|, |13[ we define a renormahzed coupling constant g^^ by the value of the four- 
point function for vanishing external momenta: 

g^R^--=-Z!Tl^'"\{Ohmo,go). (14) 
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Together with ( |T2|) one can invert the relations between m^, g'^^ on the one hand and 
mo, go on the other hand so that the (in D = 3) dimensionless renormahzed couphng 

(4) 

un := ^ (15) 
becomes the natural expansion variable of our perturbation series. Moreover with 

^:=--r?'°H{0};mo,^7o) (16) 
and ([l^) we define for later purposes 



^0 Zi{ur) 

where it is indicated that the dimensionless coupling u is to be read as a function of ur. 
To summarize, the renormalization is fixed by the conditions 



r^'°''(0;mij,Mij) = -m^ (18a) 

= -1 (18b) 

=0 

r^R'"\{0};mR,UR) = -mt^'uR = -g^i'^ (18c) 

rg'^)({0;0};mK,M^) = -1 (18d) 

(18e) 

for the renormalized (n, /)-vertex functions 

T'^^''\{p-q}-mR,UR) = [Z,{uRp-\Z,{uR)fT^:''\{p-q}-mo,go). (19) 

Broken symmetry phase 



In the phase of broken symmetry we follow |14, Q and define the renormalized coupling 



constant gR by the vacuum expectation value v of the field (p in (|^). If G^J-'^^ stands for 
the non-vanishing one-point function of the field 0o in (^1) one has 

v = vo + G'^^''\ vo = ^JSml/go ; vr:=^v (20) 

V ^3 

with a renormalized vacuum expectation value vr. We set 

^?h:=-/, (21) 
and form a (in D = 3) dimensionless renormalized coupling according to 



Proceeding as in the previous subsection we find with 

Z4: 



rriR 



1 + 



^0 



3m?, ^° , 



(23) 



and (|2ll) an expression for the dimensionless quantity 



go Zl{up.) 

u-= = ur— — - 

niR Z3{ur) 



(24) 



in terms of renormahzation constants, which again are assumed to be (analytic) functions 
of ur alone. 

The renormahzation conditions fixing this scheme are: 



r2'°H0;m^,M/j) 
Tr\p;mR,UR) 



-m 



9 r(2'0) 

Qp2 



p2=0 

3m| 



R 



i-D 

rriR Ur = Qr 



rg'^H{0;0};m«,n^) 



-1 . 



(25a) 
(25b) 

(25c) 
(25d) 



3 Determination of f+/f- 

Now we derive an expression for the amplitude ratio f+/f- in terms of analytic and 
dimensionless functions calculable by renormalized perturbation theory. 

First of all we have to realize the temperature dependence in the parameters of this 
theory. At first sight the occurrence of spontaneous symmetry breaking in the Lagrangian 
(|^) suggests t oc mQ, where according to (0) and (^ 

mo+ = ml (26a) 
ml- = -2ml. (26b) 

But when including the perturbative corrections one observes a mass shift by an amount 
ml^ in the temperature variable, which is perturbatively not calculable [jll], |12|. (This 
corresponds to the fact that for D < 4 mean field theory is no longer valid.) Distinguishing 
the different phases we have 



mR+ 



e+ ~ Ut-\ U = « - ml^) ^^^^ > (27a 



— = ^ ~ /-^:^ t- = -{ml- ml) I > (27b) 



mR 

with = t for t > and t_ = —t for t < 0. 
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In order to circumvent this problem we eliminate the bare mass mo by introducing 
the functions 

(28) 



F,(.,.):-*"* 



dml± 



go 



which are partial derivatives with respect to m"^ at fixed Qq and depend in renormalized 
perturbation theory on m/j± only. With the identity 



d_ 

dt 



d 



90 



dml 



90 



one finds 



dm\+ 




dmj^+ 






dt+ 


90 


dml 


90 




dm\- 




dm\. 






9t„ 


90 


dml 




90 



dm\+ 
dml+ 



90 



2 9^1- 



dm^_ 



(29) 

(30a) 
(30b) 



90 



Then the differentiation of ( |27a| ) and ( |27b| ) yields 

^F.{uR-) dmlJdt^ 



F+{ur+) dm\+/dt 



90 



J- 



2u-l 



(31) 



On the other hand the definition of the correlation length as the inverse renormalized 
mass implies 



f± 

J-. 



rriR- 
mn^ 



and one concludes 



F+{UR+) 



rriR- 
rriR^ 



(32) 
(33) 



So far t+ and t^, as well as mR+ and mR-, are independent of each other. In order to 
specify the approach to the critical point from both sides we choose pairs of points (t+, t_) 
in the phase diagram such that 

mR+ = mR- . (34) 
This leads to a (yet unknown) dependence between ur+ and Ur-. Consequently (|27aj) , 

F-{ur-) 



(|27bD and (|3|) read 



F+{UR+] 



and we combine these equations into the formula 

, F-{UR- 



/- 



F+{ur^ 



(35) 



(36) 



The crucial point is now that one can express (|36D as a function of a single dimensionless 
renormalized coupling constant ur, which is related to Ur± in a definite way. 
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As mentioned earlier, the critical theory is characterized by the non-trivial fixed point 
m|j of ur. This is equal to the non-vanishing zero of the renormalization group /5-function 



P{ur) := niR 



d 



dm 



R 



ur = - 



d 



90 



du 



R 



-1 



ln(M) 



D = 3, 



(37) 



according to m = u{ur) = go/rriR. Thus, if we are looking for a new coupling ur, which is 
to be used as a renormalized coupling in both phases, we have to ensure that this coupling 
leads to the same /3-function in both phases. Consequently the functional relation between 
mR+ and ur in the symmetric phase must be the same as the one between rriR- and ur 



in the phase of broken symmetry. In view of (|37|) one can comply with this condition by 
using the relations (|T^ and (0), which represent perturbative expansions of u± in terms 
of Ur±: 

U- = h^iuR-) (38a) 
u+ = h+{uR+) . (38b) 

One possibility is to identify ur with ur+ in the symmetric phase and to define it in the 
phase with broken symmetry by means of 

u_ = h+{uR) . (39) 

The second possibility considered in this work is to identify ur with ur- in the phase 
with broken symmetry and to define it in the symmetric phase by means of 

u+ = h^UR) . (40) 

In both cases the /3-functions relevant for the high- and low-temperature phases 

d 



P±{ur) 



duR 



In ('u±( 



Ur 



are equal to each other, 

13+{ur) = p-iuR) 
In the first case their common fixed point value is 



u 



R+ ) 



for the second choice it is 



With the definition 



u 



R 



U 



R- 



^{ur) 



-\Ur~ {Ur, 



Ur+ [Ur 



(41) 
(42) 

(43) 
(44) 

(45) 



[2$(u 



R) 



(46) 



we can evaluate ( p6D at the critical point, 

l± 

f- 

to get numerical values for the amplitude ratio /+ / /_ up to a given order in perturbation 
theory. 
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4 Results 



We have calculated the renormalization functions necessary for the procedure outlined 
in the previous section in perturbation theory up to the order of two loops. In the 
symmetric phase 7 massive Feynman graphs and in the phase with broken symmetry 24 
additional graphs have been evaluated analytically in D = 3 dimensions with dimensional 
regularization. We omit the details of our calculations and present only the main results. 
In the symmetric phase we obtain the bare expansions 



m 



R+ — 



(4) 

9r = 



0+ 



Ur+ 



90 



1 - ^ Mo+ ^ 
Stt 

3 



1 



167r 
1 



647r2 
5 

162^ 



Ur 



79 
162 



3 



(div) 



ul+ + 0{ul+)^ 



216 6 



(div) 



+ 0{u, 



0+ 



'} 



0+, 



(47a) 
(47b) 

(47c) 



which are inverted to 



ml+ 



Ur+ 



139_1 

216 6 ^+ 



245 
[162 



is 



(div) 
1+ 

(div) 



R+ 



U+ = Ur+ < 1 + 



Ur+ + 



575 



167r 207367r2 
In the broken symmetry phase one has 



4+ + 0(4+)} 



+ 0(4+ 
+ o(4+)} 



(48a) 
(48b) 
(48c) 



gR 

ur- 



m. 



0- 



1 



3 



647r 



90 

Uo 



7 
31 



Uo- + 



1287r 



1 2 
+ 647r2 ""0^ 

37835 
331776^2 ' 

1 o 



19525 
5184 



3 



(div) 



647r2 



u 



0- 



m25 



13824 



3 



(div) 



+ 0{ul) 



+ 



0{ul)] 



(49a) 
(49b) 
(49c) 



and 



uo- 



U- 



ur- 1 1 



31 



1 



Ur- 



nin- — m 



R- 



1287r 647r 

3 1 
Ur- - 



Ur- 



647r 



2 



647r: 
2191 

1658887r2 



23663 1 „ 
13824 3 ^ 



45125 2 

h --B 

10368 3 



R- 

(div) 
R- 



U 



R- 



+ 0{u%-) 



+ 0(4-)} 
+ o(4_) 



(50a) 
(50b) 
(50c) 



The divergent pieces 



ry(div) 1 , (ml±\ , const. 

si --InU- -7+^ + 0(e) (51a) 



d(<^''') 1 1 f^R±\ , const. 

7-1^^ -7+^^ + 0(6) (51b) 
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vanish after differentiation with respect to mQ±. So we get from (|28|) and the above 
equations the finite functions 

= -1-^4^ +0(4+) (52) 
F_ (ur- ) = 1 + ur- =^ ul- + O iul- ). (53) 

As discussed in the previous section the new couphng constant ur is to be introduced 
now. We consider both choices explained above. 

1. Expansion in Ur+: 

We determine ur~ as a function of ur with the help of equation (^) . From ( [48c|) 
and (|50c| ) the result is 



/ \ 9059 \ 

= (l - 32^ + I6588S5 + ' 

Since m|j = u*^+ we identify m^j with ur+. Then 

<I>+(mr+) = 1 + H ^^-^4+ + 0(4+) • (55) 

+ ^ " ^ 1287r ^ 163847r2 ^+ ^ ^ ^ 

2. Expansion in Mr-: 

We adjust ur to the low temperature coupling Ur-. The equations (|40|), ([48c|) and 
( ^Oc|) show that in this case depends on ur via 

= Ur(i -\ Ur ^ Mr + 0(Ur)^ , (56) 

^ " V 327r " 1658887r2 « ^ ^ ' 

and = allows the identification ur = ur-- Consequently 

^(u^-) = 1^-^ur-^ ^^-^4- + 0(4-) • (57) 

^ " ' 1287r " 163847r2 ^ ^ ^ ' ^ ' 



In order to get numerical results for the amplitude ratio f+/f~ according to 
we need the values of and u. The high temperature fixed point and the index 
z/ have been calculated in the framework of renormalized perturbation theory in D = 3 
dimensions in [|1^, . Therefore the whole analysis only requires the present methods. To 
improve the quality of the results we prefer, however, to make use of additional information 
available about and u. For the exponent i' there are recent Monte Carlo result in 
| T6| . The fixed points have been estimated from the longest presently available high- 
and low-temperature series in |T^, and for we also use an estimate cited in 0. To 



estimate the sensitivity of /+//- to these parameters we list the values resulting from 
different choices in tables 1 and 2. 

The functions $±(mr±), given in equations (|55D and (|57D , have been evaluated at the 
fixed points with Fade approximants. By this we do not mean that the perturbative 
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Input 


Amplitude ratio f+/ f- 


Exponent v 


Fixed point 'U^+ 


[2,0]-Fade 


[l,l]-Fade 


[0,2]-Fade 


0.6300(15)0 


23.73(8)^ 


2.2348 


2.2663 


2.0716 


24.56(10)0 


2.2605 


2.2956 


2.0775 


0.624(2)0 


23.73(8) 


2.2178 


2.2487 


2.0570 


24.56(10) 


2.2431 


2.2775 


2.0630 



Table 1: Ratio /+//- for given values of v and 



Input 


Amplitude ratio f+l f- 


Exponent v 


Fixed point m^.. 


[2,0]-Fade 


[l,l]-Fade 


[0,2]-Fade 


0.6300(15) 


14.73(14)llZl 


2.0122 


2.0497 


2.0329 


15.1(1.3)1 


2.0255 


2.0660 


2.0546 


0.624(2) 


14.73(14) 


1.9988 


2.0356 


2.0255 


15.1(1.3) 


2.0119 


2.0517 


2.0406 



Table 2: Ratio /+//- for given values of v and 



expansion is analytic, which it is presumably not, but just take the different results as a 
measure for the error due to the shortness of the series. The application of Fade-Borel 
methods, which we also tried, does not yield any improvement. 

We summarize the contents of the tables by extracting mean values of /+ / /_ for the 
high and low temperature fixed points: 



2.18(12), U± =2.03(4). (5^ 





The number given in brackets is the maximal variation owing to the inputs of u*^± , v and 
the different Fade approximants. 

The fact that the variation between the different Fade approximations in table 2 is 
much smaller than in table 1 indicates that the low temperature renormalized coupling 
constant ur- is the better expansion variable for our series. An examination of (p3| ) and 
(|5^) confirms this impression since the first and second order perturbative corrections to 
the leading terms in the coupling ur- are smaller than those in the high temperature 
coupling uji+ (41% and 8% against 66% and 9%). Therefore we tend to have more faith 
in the estimate resulting from the expansion in ur- . 

For the sake of completeness we add a remark about the perturbative determination 
of the critical exponent v. It is related to the renormalization group functions r/j, for 
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i = 1,2 defined by 

d 



90 



\n[z,±{un±)), (59) 



through the equations |]T| 
1 



2~7]3±{ur±) +r]2±iuR±) , iy = i^±iuR±). (60) 



The renormalization constants Zi, i = 1,2, are calculated straightforwardly from (|47a|) - 
( |50c| ). The expansions of which are of the form 

i^± = l + 0{uR±), (61) 

can be inserted into (^) and evaluated at the fixed points u*j^± in the same way as before. 
The expansion has the form f+/ f- = V2 + 0{ur±). This is the typical structure for 
renormalized perturbation theory because the higher order corrections in ur± represent 
the in D < 4 non- negligible deviations from the known mean field value (/+//_ )m/ = \/2. 

With this method we find {f+/f^)ht = 2.22(6) and {f+/ f-)u = 1.86(6). We must 
emphasize, however, that these results are strongly influenced by the behaviour of the 
short series for i/± and therefore are certainly less reliable than the results above, where 
more precise information about u has been employed. 



5 Conclusion 

Our results for the universal ratio of correlation length amplitudes indicate a value of 

/+//_ = 2.04(4) . (62) 
The theoretical estimates in the literature are 1.91 from the e-expansion ^, and 1.96(1) 



T8| , 1.94(3) [jl^ from high- and low-temperature expansions. They are deflnitely below 
our number. On the other hand, experimental values from binary fluids are 2.05(22), 
2.22(5) [|T9l and 1.9(2), 2.0(4) [gg. 

In conclusion it appears quite reasonable to us that the correct value of the universal 
amplitude ratio of correlation lengths can lie above 2. The calculation of further or- 
ders and a detailed investigation of the resulting (non convergent, but hkely asymptotic) 
perturbation series may lead to more accurate estimates. 

Finally let us mention that with an appropriate modiflcation our methods are also 
qualifled to determine the amplitude ratios {f+/f-)ph of the 'true' correlation length 
defined via the physical mass, and C^/C^ of the susceptibility. We have performed the 
calculation of C+/C_ in the two-loop approximation, too. The result is consistent with the 
estimate in [jlT|; their series from Feynman graphs in the symmetric phase is, however, 
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longer. The calculation of the 'true' correlation length in the one-loop approximation 
reveals that the amplitude ratio is only negligibly different from the one considered here. 

As a consistency check we have also reproduced the e-expansions of the quantities 
considered in this work. 
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